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abstract 


,oltttioJM  in  graphical  form  are  pre.ented  for  the  .hock 
•  l/a  °f  a  ?ne"^e8r«e-of-freedom  .y.tem  subjected  to  some  ideal- 
devll(on4Cln^UnCt|i0"‘*  E1*m«lM*IT  procedure.  are  outlined  for 


1.  INTRODUCTION 

olacemenTo^r^"8  °f  en«jneerjn8*  knowledge  of  the  maximum  dis- 

5f, charac^Li.tiiTe111*"1  8y8tem  i8  de#ired  aa  *  function 

a  .hock Sy,tem<  Thi«  upuctrum  i.  called 
a  anocK  spectrum  (ref  1,  2)  when  transient-type  forces,  such  as  in 
an  explosion,  are  applied  to  the  system. 

*y,tem  °f  interest  in  this  report  is  a  one-deoref-of 
(ODOF)  undamped  system  with  a  linear  restoring  force,  !  e  a 
fe«nteHma*“?rln?  8ystem-  The  interest  lies  fn  the  effe’ct'of  dif- 

mtU  ly  atTeg.t  -Gr°on*H  The,8tudV  i8  limited  to  a  system 

initially  at  rest.  Ground  accelerations  and  velocities  are  not  consid- 

the  ^  irr*levant  whether  the  maximum  displacement  occurs  in 

the  P°»itive  or  the  negative  direction.  A  structural  member  eener 

In^theVituation  chos^Ue8th^  Y°Un*'8  modulu8  in  tension  and  compression. 

the  exceXs  wil  be  ;ofed.P°Slt,Ve  ™XimUm  18  U8ually  Ur«er- 


Two  chief  types  of  forcing  functions  are  considered: 

(1*  andalf"8ine  PUl*e  With  suPerimP°8ed  Dirac  impulses 
(2)  A  tJ  exp(-Pt)  function. 

Pulse  (1)  will  have  impulses  of  two  types: 


(la)  five  equispaced,  equistrength,  positive  Dirac 
impulses  with  the  first  at  time  zero  and 

(lb)  the  preceding  with  five  additional  negative  impulses 
symmetrically  interspersed  (alternating  impulses). 

A  Dirac  impulse  of  strength  F,  will  be  one  for  which 

/  F,  5(t)dt  =  Fi 


The  scope  of  this  report  is  limited  to  presenting  the  mathe¬ 
matical  details  used  in  evaluating  particular  shock  spectra.  Their 
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"h‘o*=v  ""“'1  bf  in  •  inter 

SKMSsaE,— 

pul«e.  (ref  4).  Moreover,  the  V-  *  often  ob*erved  on  .nod 

jectile,  .triicin,  .,«, 

^r,tr7"^ 


^-g£N£M.L  DEVELOPMFmt  ^  TrTTnn 
We  are  concerned  with  ODOF  system.  4 

the  equation  yetems,  i.e.,  systems  described  by 

mx"  +  kx  =  F(t) 

The  rest  conditions  assumed  are  (1) 

x(°)=0,  x'(0)  =  0 

™e  forcing  function.  r,„  con.*..,., 

n..riou..pVcT»"iV.“"!”<lT“'1',l'ir,',‘^;  °“r  ODOf  -ynten,  i.  „„ 

f*  V*n‘".  f‘8ure  2.  x  be  definedbvvt  ^  K^eralizaUon 

in  relation  -f.tbehrr:uynd^ne^reectiie^nyge8 

motion  of  the  system  may  be  written  P08mon8‘  The"  the  equation  of 


mx'  +  kx  =  F(t)  .  mz" 


IS™.'™'  13  thu*  e^'valent  to  a  force  F(t,  of  opposite  sign' 


(3) 


Figure  i. 


Physical  model  of 


generalized  ODOF 


system. 
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Normalization  of  the  variable*  in  (1)  will  have  two  advantage*: 

(1)  the  equation*  to  be  *olved  are  easier  to  manipulate,  having  a 
minimum  number  of  e**ential  variable*  and  parameters;  (2)  the 
solution*  represent  many  cases  resulting  in  economy  of  computation 
and  representation.  The  first  step  is  to  normalize  these  variables. 

One  may  define 

F(t)  =Fuf(t)  (4) 

where  F0  is  the  peak  value  of  the  forcing  function  excluding  impulse 
components  (fig.  1).  V  we  further  defire 

,:N?  s  k/m»  xs  =  Fo/k,  u  =  x/x,  (5) 

equation  (1)  reduce*  to 

u"  /'»n2  +  u  =  f(t)  (6) 

Derivatives  so  far  have  been  with  respect  to  l.  In  (6),  nj^j  is  the 
natural  frequency  of  the  system,  xg  the  static  displacement,  and  u  the 
response  factor.  We  have  yet  to  normalize  ti  ■ne;  one  may,  in  general, 
define  a  unit  of  time  tc  and  let 

r  =  t/tc  (?) 

Examining  figures  la  and  lb,  we  see  that  it  is  most  natural  to  define 

(8) 


':Flc 


where  i  p  is  a  frequency  characteristic  of  the  driving  function.  To 
allow  (8)  to  represent  the  situation  in  figure  lc,  we  may  allow 

0  =  up/w  =  l/tc  (9) 

The  characteristic  forcing  frequency  iB  thus  x@  and  th?  peak  value 
occurs  at  2tc.  Furthermore,  we  normalize  the  natural  frequency  as 

V  =  Wjq/'l’F  ( I®) 

With  these  relations,  (b)  becomes 

u '  (T)/vP'f:’  +  u(T)  =  f(-r)  (11) 

Without  ambiguity,  derivatives  of  normalized  variables  are  with  respect 
to  t.  In  addition,  some  liberty  has  been  taken  with  the  symbols  u  and 
f;  as  functions,  they  are  not  the  same  functions  of  t  that  they  are  of  t. 
The  symbols  in  (11)  will,  with  no  real  ambiguity,  represent  the  new 
functions.  See  appendix  D  for  a  table  of  symbols. 

The  next  step  in  the  analysis  is  the  solution  of  (11)  for  u.  For 
the  functions  presented  in  figure  1,  the  solution  may  be  carried  out 
explicitly.  There  are  several  known  ways  of  achieving  this.  The 
author  has  chosen  the  Laplace  transform  method,  but  any  of  the  other 
standard  methods  are  equa  Uy  quick. 
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of  the  norm^ h^Lmna^Ta7freqSuPency  T"  For  I,'1' f°r  each  val“« 

'•  the  exphcit  8oiution  u  may^be  aj^e„^er„r"vrdpieB  function‘ 

UM  =  u(-M)  M 

cases  because  ^of  th  ^present-1  (  f  =  °i*  ThU  *"  not  feasible  in  our 
to  be  no  better  recourse  than  to  .cTn  l'1!!'  ?ea*8  in  u<  There  appears 
of  values  of  t,  and  find  the  maximax  fthl*81!  kT  f°r  4  8uilable  number 

comparison.  ^  in  refeVencrV^a  i^ne8  ncsative^ens^^by  emUm^ 

with  experimental  forcing  functions"  ThPr°hram  **  described  to  deal 
program  was  not  easy  to  make  bu?  the  ch°lce  towrite  a"  ad  hoc 
aT.ti°n#.analytically.  the  general  c^um??6  °f  fl?din8  the  resP°nse 
it".  favorV.,rtUe  °f  added  in81«ht  w*th  the  a dlToc  program  Sir8"8"1’’ 
Except  for  large  values  of  v  , 

jr'~  s  szzzzjz:  ich 

not  seem  worthwhile  to  improve  the  accTra^  °h  th<?  results-  il  does 
gence  algorithms.  F  accuracy  by  use  of  local  conver- 

*"  •»  (Ming  „M  ,h. 

free-vibration  era  (residual  9rJ?r  y.*  also>  to  find  the  u,J1n  the 
-ay  be  larger  than the ^inZYoTL**  f°r  CTerta’"  val-s  of  it 
ransient-like  functions  (cf  fig  lc\  no  f  era* .  *n  tbe  case  ot  so—e 

strictly  speaking.  I„  tho.e  1“  :-vlbration  era  exists, 

-ay  be  simply  derived  from  the  forced  d°eS  exist’  the  “Pectrum 
damping  is  absent.  Thus  Trom  in a  Cra  resPon8e  factor  u.  when 
we  must  have  ’  consideration  of  constant  energy, 

7kx?(tc)  +jmx\)  =  7  kx  wr 

where  t  is  chosen  to  be  t  u  *  M  (12) 

to  the  system  by  the  external  forcT  The  T*  ent'rgy  haS  been  added 
to  the  maximum  potential  enercv  in’. hr  energy  is  then  equated 

nation  of  (12,  results  in  the  expression"  free" vibral-"  era.  NoLalT- 

. .  ,  .  um  '  "^u -  ( i )  +  ru'(i,/YTr,c 

:^rssr-t8-y.  uM(,J) 

of  obtainmg  the  free  era  spectrum  by  roans^f’  th 'V ltefna,iv“  m&od 

be  mnl *“,  related  by  a  constant  factor™ the  fi  °UrK'r  "P^trum 
be  more  direct.  lo  th‘  free  spectrum  _  may 

sr^ss:sssL^ 
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3.  HALF-SINE-PULSE-TYPE  FORCING  FUNCTIONS 

Aa  a  firat  specialization  of  equation  (11),  let  ua  consider  driving 
functions  of  the  type  given  in  figures  la  and  lb.  Since  we  are  dealing 
with  a  linear  system,  we  may  find  the  responses  to  the  pure  half- sine 
pulse  and  the  delta  function  train  separately,  then  add  the  two  pro¬ 
portionately  to  obtain  the  total  response.  (We  may  not,  however,  add 
their  separate  spectra  to  find  the  total  spectrum. ) 


Since  we  use  (13)  to  obtain  u^  directly  in  the  free  era,  the  response 
is  needed  only  for  t  s  1.  The  normalized  half-sine  pulse  may  therefore 
be  written  as 

fHs(T)  =  sin  tit  (1^) 

and  the  equation  of  motion  reads 

u"/v2w2+  u  =  sin  tit  (1®) 

By  standard  methods  we  obtain  from  (15)  the  transfer  function 

UHS<«>  88  (.8  +fi(.'WiT7T  (16) 


UHS'*‘  -  («“  +1PW7TW) 

With  the  aid  of  a  partial-fraction  expansion,  one  obtains  the  desired 
inverse  transform 


uHs<t)  =  I'-jV.?  (sin  VTW  •  V*inTir) 


Next,  for  the  finite  train  of  positive  delta  functions,  we  have 

fo+(T)  =  MT)  +  MT -l4)  +  ’  •  •+  5(f  -(m',*4)  (18> 

whose  transform  is 

Fd+(s)  *  !♦•-•*  +  •  *  •  +  */s  (19) 

The  function  fD+  represents  the  situation  after  the  m,h  impulse  and 
before  the  m+nh.  where  m  =  1,  2,  3,  4,  5.  The  train  of  negative 
impulses  is  represented  by 

fD-<T>  =  -WT-  V«'>  (20> 

and  has  a  transform 

FD»  =  -e'8/,n  Fj^s)  <21> 

It  follows  that  the  transfer  functions  for  these  cases  — the  left  side 
of  (15)  being  still  valid  —  are 

U^s)  =  fv=V/(s2  +  Y2tr2)}Fm(.)  <22> 

UDJ.|.  »wl-l  |i3> 

whose  inverse  transforms  become 

ud^1t)  =  -vtf®ittVtTr+  sinv(T  -  Vb)tt  +  . . .  +  sin  v(T*  (m+  "  (2^ 
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And 

“DJT’  -  . in v(T-(2m. 

m+e=' mT+or,3.”*a^?I^  °‘  ”  *"  (19)  and  al«°  either 

bined^  ^ZlLTtri^^  in  <2<>  ~  (25)  may  be  com- 

However,  IT'  J  >*.«<%> 

never  greater  than  five,  and  the  °£  terms  in  the  Bum  is 

pole,  for  ceruin  value,  of  Y  which  ^ecui"™  artlficiallV  introduce, 
digital  program.  Y  require  .eparate  handling  in  a 

•in*  JSu.’JiS.-  “»  '••**<!  **.*,„.. 

u  =  UHS  +  r(«d+  +  Lud_) 

amplitude  of  the  pur e^alf. Vine' pull*"' Lh-°0  f  function  to  the 

-  ,h“  »•  vs^asrzssssss.^ 

“■  "“‘‘""’n-nnun  SJ 

“M  *  Rr"w"’  *  ,,3*, 

.  “ub,tltut,°n  of  (17).  1241.  /psi  . 


the  substitution  of  (17)  1241  /7ei 

^ir!^0  (13a)  oi't  In  ^ated 

—  FORciNG  FUNCTION  of  thf.  tvpe  taex^-ft^ 

The  t-dependent  function  may  be  written  r  k 

F<„.C*,V».  Qj  being  an  arbitrary  constant,  as 

The  peak  value  of  F(t)  i.  found  to  be  (27) 

F:  =  4C0/@V 

Consequently  (28) 

and  f(t)=i^t8e-6t 

f(  t  )  =  ^-tVt  29> 

when  we  use  (7,  and  (9).  which  together  imply  that  t  .  8t.  ^ 
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It  i*  con- 


The  usual  liberties  have  been  taken  with  the  symbol  £. 
venient  to  rescale  u  and  y  as 

u  =  eaw /4,  ?  =  yw 

Then  equation  (11)  becomes 

w"/?2  +  w  =  T3e“T 

Using  the  Laplace  transform,  one  readily  obtains 

=  {s^ggj('s~TF 

One  may  proceed  now  by  two  different  paths.  The  partial  fraction 
expansion  may  be  effected  and  the  inverse  transform  obtained  from 
the  fragments;  or  the  convolution  theorem  may  be  applied  followe 
by  an  integration  by  parts. 


(30) 

(31) 

(32) 


4.  1  Partial  Fraction  Method 


Writing 


7^+7Tl+(7TTF  +  (ITI? 


(sg+?g)(s  +  l? 

one  obtains  by  standard  methods  (apx  A) 


.  e2-  3 

A"(l  + 

_  l  -  3e;g 

B  =ttt?!F 


(33) 

(34a) 

(34b) 


and 


C  =  -A 

D  =  2/(1 +  f)3 
E  =  1  /( 1  +  ?z) 


The  solution  follows  by  substitution  of  these  coefficients  into  the  in¬ 
verse  transform  obtained  from  (32)  and  (33): 


(34c) 

(34d) 

(34e) 


w(t)  =  2f3(Aco8  ft  +  (B/?)  sin  ?t  +  (C  +  D-r  +E2T2/2)e*T} 


(35) 


4.  2  Convolution  Method 

Applying  the  convolution  theorem  to  equation  (32), 
solution  to  (31)  in  the  form 
T 

w(t)  =  5 fvae“vsin  ff (t  -  v)dv 
0 


one  obtains  the 


(36) 
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Mor,“%«',r,u"°"„Th'“«  ?„"*  ";y/°nv*«  n‘K7,  ta™  o?T35)on  lh' 

For  reference,  the  explicit  solution  is  given  below  as 
W<T,=n+^F{C<I-3?a)»in?T-5(3-?a)co.5T]+5[(3.5=)  +  2{1+52)T  + 

!4(l+(r®)3T2]e*T }  (37) 

—  PROGRAM  FOR  half-sine- pulse  functions 

The  FORTRAN  IV  program  does  the  following: 

tabled.  (1>  ThC  resP°nse  fact°t  u  is  calculated  for  0  <  T  <  1  and 

determine  SSSB?,!?,. 

ized  natural  frequency!  UM'S  t3bled  38  3  function  Y.  the  normal. 

(4)  Tne  positive  u^s  are  graphed  against  y. 

(5)  The  residual  uM>s  are  tabled  and  graphed. 

which,  among  other  posTibilfues Utallow°th8  incorporated 

selected  individually.  These  are’  note^T-*^  precedlnS  items  to  be 
program  MAXU  in  figure  3  The«  d  t>rlefly  ln  the  comments  to  the 
the  positive  impulse  case  c'an  be  s!lec°tm^Knt,8  8h0W'  mo-°ver,  that 
that  special  test  value,  of  v  mav  *ettin«  M1  be  "°^nity; 

M2  =  1);  and  that  tables  for  a  redded  n  M5  =  1  <U8ed  with 

chosen  when  KTEST  <  KK  The  Mrt’  n"mber  of  R  values  may  be 
represented  by  the  KT  a/ray!  partlcular  of  R  desired  are 

contairdirrheH^  s1b7ouUnempZTOpPi0ttteH  **  U'in«  the  featu~ 

is  read  in;  for  each  plot  (indexed  bv  Kl^uT'  ,the  number  of  plots 
ICURVE  array;  these  values  represen  'th^  M  “7  ^  read  into  the 
make  up  a  family  of  curves.  pre8ent  the  M  values  of  R  chosen  to 

scales  UWar^ly  of curveTfoTetch t/aph!  A-.Hausner's  SCALE2*  — 

iCnd7pee8ndentnoef  Ihe'five?6"  rC'Pl0t  thl'6e  °f  ^^^whfclTaL15^^ 


"^8eVer3i  8ubrou‘>n-  documented  at  HDL  for  internal 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


DIMENSION  TAUItCOOl.TAUPUlOOOI.STAUPIllOOOl.SGTPKlOOOl.UI  1000.51 
l,CTAUPII1000liUNAXISI.UM|NISI.R!5).KT(5I.URI5l  iCCHIO I .**!»•» 

2  UMINAISI.UMAXAIS! 

OINENSION  GtSOOl.UNISOO.SI.URESNISCO.SI 
i  FORMAT!  14 19 ) 

«  FORMAT!  20X4MR  ■  F7«4ilOXIH61NNA  ■  Fo.S.lOXTHUMAX  ■  IPElS.Ti 
1  10 X7HUMIN  ■  1PE1S.7//10X4(TX1MT12X1MU9XII 
12  FORMAT! II0X4I4X0PFA.S.IPE29. 711 1 
22  FORMAT!  12X5NGANNA5I 7X7HUR4X  R.F4.4II 
J2  F0AHAT(/nx0PFA.2.1P5C29.TI 

J2  FORMATI^SOXJRMPOSITIVE  IMPULSES  AESIOUAL  SPECTRUM10X4MR  •  , 

1  F7.4///I 

42  FORMAT!  I1SXSI2X5HGANMAIX5MURMAX4XI I  » 

72  FORMAT! (15X3(9PF00.1.1PE14. Till 
42  FOAAATI2OX25MFORCE0  VIBRATION  SPSCTRUN//I 

*2  FORMAT!  2CX4SHTEN  ALTERNATING  IMPULSES  ON  A  M4LF-SINE  PULSE//) 

102  FORMAT! 20 X4SNFIVE  POSITIVE  IMPULSES  CN  A  HALF-SINE  PULSE//I 

112  FORMAT!  J9X42H  ALTERNATING  IMPULSES  RESICUAL  SPECTRUM10X4HR  •  • 

1  PT.4///I 

122  FORMATIITXSITxTmUMIN  R.F4.4II 
192  FORMAT!  IF II 


II 

JJ 

RR 

Ml 

P2 

MS 

M4 

MS 


M7 
RTEST 


nlmoer  of  tau  values 
nlmber  of  gamma  VALUES 
nlmoer  of  r  VALUES 

ALTERNATING  IMPULSES 

NRITE  TABLE  OF  U  VS  TAU  IFORCEO  ERA  I 

VRITE  TABLE  OF  UMAXtUMIN  VS  GAMMA  IFORCEO  ERAI 

PLOT  UMAX  VS  GAMMA  IFORCEO  ERA) 

READ  GAMMA  IN  IREAD  CARO 

NRITE  TABLE  OF  URESM  VS  GAMMA  IRESIOUAl  ERAI 
PLOT  URESM  VS  GAMMA  (RESICUAL  ERAI 
.  NUMBER  OF  VALUES  OF  R  FOR  nMICH  TABLES  ARE  WANTED 


REAOISf 1 1  II.JJ.RR  .  .. 

READ! Si  II  M1,M2,M3,M4,PS.M6iM7,RTEST 
REAOIStll I  OTAUiOGi (Rl* I iR*l «RR I 
AEA0IS.1I  IKT1KIiR«1iKTESTI 
PI  •  1.1415427 
EL  •  Ml 

IF(Ml.NE.l)  EL*C. 

DO  119  RA-l.RTEST 
R  •  RTIRA I 
119  RAIRAI  •  R(R) 

CO  9  1*1.11 

FI  •  1 

TAUIII  •  FI.CTAU 
TAUPIIII  •  PI.TAUIII 
STAUP1 1 1 )  ■  SIKTAUPIIIII 
(TAUPIIII  •  COSITAUPKII) 

9  CONTINUE 

IFIMS.EO.il  CO  TO  40 
00  19  J'liJJ 


Figure  3.  Program  MAXU  and  subroutine  PL.TOPT. 
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PJ  m  J 

6(JI  •  FJ«OC 

i*  continue 

CO  TO  TO 

‘  „  "*'*■“•'  ‘“"Hit  «  „ 

S  S'K’-ili.r 

unm*ki  *J*o?o 

UNINIKI  .  0  « 

>*  CONTINUE 

00  •*  M,10 

PH  •  M 

JCP  ■  FN. Cl/10. 

521?*  •  C0$(F5P| 

•*  CONTINUE 
00  **  Mill 

5SZri",.;.s,N(c‘j»*T*u^niM 

mi.oI.iooTm’tS 

MICK  .  0.0 
CO  TO  1C 

15  I?f  *;6T »*00  J  GO  TU  25 

co  7o*fo*“  00  T0  110 
110  cT5S  ;0«®"*tt-i»3i 

**  JH  I.GT.300)  GO  TO  35 

loin  It'"  00  10120 

120  ao'ro  ;oSCT>,,'-»oo» 

”  mi;  “I;4?”  00  T0  45 
go 7S*Io  50  T0  130 

l"S,S.VW,WM»  *  5GTPI |  1-300) 

45 1 £» cIa5?0*  Cf1  T°  « 
is  ?i*3o  60  T0  140 

140c#S,?S3*o$6m',-,M»4«"m.,3oi 

**  ml:5!;4?01  50  TO  65 

CO*TO*lj’*  C°  T°  l,° 

*"  S'S 

45  /rjirGT •,03 *  00  TO  T5 
CO  TO**?***  G°  T°  160 

S'S 

Figure  3,  Program  MAXlt  =  . 

and  subroutine  PLTOPT. 
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IFINi.EQ.il  CO  TO  1T9 
110  MICK  t°SGTPlll-100l  ♦  SGTPI I 1-3901  ♦  SGTPttI-5001  ♦  SGTPIII-TOOI 
GO  TO  40 

•)  IFtl.GT.990l  GO  TO  99 
IFtNl.EQ.il  GO  TO  180 

180  BRICK  ■  SgTPIII-IOOI  ♦  SGTPIIt-3001  ♦  SGTFItl-SOOl  ♦  SGTPItl-TOOl 
GO  TO  50 

95  IFINl.EQ.il  GO  TU  190 

190  «ic2  «°$GTF1 I t-1031  ♦  SGTFltl-3901  ♦  SGTPIII-SG01  ♦  SGTFII1-TOOI 
1  ♦  SGTPI 1 1-9901 
GO  TO  50 

10  BRACK  •  SGTPI  (I  I 
GO  TO  103 

20  BRACK  •  SGTPI til  ♦  SGTPI 1 1-200 1 

)0  BRACK  i°SGTPItII  ♦  SGTPltl-2091  ♦  SGTPIII-4001 

40  BRACK  i°SGTPI(ll  ♦  SGTFItI-2001  ♦  SGTPIII-4001  ♦  SGTPItI-6001 

50  MA«  »°SGTP  III)  4  SGTPI  1 1-2001  ♦  SGTPI  tI-4001  ♦  SGTPI  tl-BOOl  ♦ 

1  SGTPI 1 1-890 1 

100  ^81  ABSt  Gt  j7-l .9 1  .LT.t.E-51  HSO.5.1  STAUPl 1 1 1-TAUPI 1 1  LCTAUPI 1 1 1 1 
DEL  TP  •  C2*BRACK 
DEL TN  ■  C2*BRICK 
00  59  K>liKK 

UlltKI  «  HS  ♦  RIKMIOEITP  -  eu.DELTNI 
IFtUI  I.KI.GE.UNAX1KI1  UNAXIKI.Utl.K) 

IFtUt I.KI.LE.UNINIKII  UNlNtKI-Utl.RI 
INIJ.KI  •  UNAXtKI 
59  CONTINUE 
49  CONTINUE 

SStJl  /cUlCGPtlOlU.I  ♦  C2.RIK»pICGPU0IaCGPIBI4CGPI6I*CGPI4I* 
1  CG5t2l-EL.tCG»(9|.CCPtTl*CGPt5|4CGPtJI«CGPIllll 
URESNIJtRI  •  SQRTtUt II,KI»Ut 1 1 ,KI  ♦  URIXl»URtKII 
99  CONTINUE 

DO  69  KA-l.KTEST 
K  •  KTIKAI 
UNAX  A  (RAJ  •  UNAX  t  X  I 
UNI NA IK Al  •  UPIMKI 
IFIN2.NE.il  GO  TO  69 
IC  *  0 

T  6RI  TE 16.132  I 

IFIMl.EQ.il  NRIT6IA.92I 
IFIHl.NE.il  NRITE  16.102 1 

Ml  tI  tt!**'  R  IK1,  GIJI.  UNAXtKI,  UN  INI  K! 

II  •  1  ♦  IC*200 

III  TE  16 .12^  I  TAUI 1 1.Ut  I  .Kl.TAUt  I*5C  I  »Ut  I*50»XI  .TAUl  I*1CCI» 

1  Ut  1*130  iK  I  ,  TAIIt  1*1501  ,Ut  1«159,KI,  I  *11, 121 


Figure  3.  Program  MAXU  and  subroutine  PL.TOPT 


c 

c 

c 


GO  TO  69 

GO  TO  T  *  * 

*«  CONTINUE 

JTIHI.NE.ll  CO  TO  19 

“  -»■  •> 

NR1T£<6,*2| 

If { 2* •*«•*»  “*mt6.92» 

ITIHl.NE.il  ufl  I  TK  (a  .  106  1 

39  ‘^®»af:fisn 

IT(H6.N£.1|  CO  TO  12J 
CO  109  HA«1.KT£ST 
*  •  KT(KA| 

ICE  •  o 

200  U  •  1  ♦  ISO* ICE 
L2  •  Li  «  49 
NR! TE 16 (1)2 ) 

!» Hi «<•<» 

NRI Tt (6 .42 1  MrTEI6'^»  *<*» 

. . 

109  cKL"*’,G0  T°  200 


'9  |  r9V»  ' 


*  two 

IT(HT.‘;S.'{J  call  pitopm 

REMIND  9  *■  opt*3»G.URESH,  JJ,kk  | 

STOP 

EM} 


1  TCRH?n?6r»£;,S"V,>,,’<’'CU,,/?'5-R'*«NI';.*..NAI« 

N  •  JJ 

REAOIS.l)  NOP T  S 
OC  9  KM.NQPlTS 
«fAOI5,n  H 

t  )  C^LL  S'  ilCA*/  g  u  4 

TTIL.EO.^1  call  srflfj  f'f’I' 

9  CONTINUE  :>  ■*ie6l*,Y,7>|o,2t2,u, 

RETURN 

END 


Figure  3.  Program  MAXU  and  subroutine 


pltopt. 


*  Gl L ♦ ICO  I , 
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The  calculations  in  the  bulk  of  the  program  are  done  so  as  to 
avoid  repetitious  calculations  of  the  sine  and  cosine;  arrays  such  as 
TAUPI,  STAUPI,  SGTP1  are  precalculated,  keeping  the  running  time 
to  a  reasonable  5  to  8  min  on  the  7094,  A  bettei  method,  though, 
uses  recurrence  formulas  (ref  8)  to  generate  all  but  the  first  members 
of  the  circular  functions  with  multiplications. 

When  y  =  1,  equation  (17)  has  a  pole;  thus  Cn  after  statement  29, 
becomes  zero  on  the  7094;  HS  in  100  is  thus  zero,  but  the  IF  state¬ 
ment  following  it  —  allowing  for  binary  roundoff  —  calculates  the 
limiting  value  obtained  by  L'Hopital's  rule. 

Many  of  the  variables  have  names  suggestive  of  their  counter¬ 
parts  in  the  previous  sections;  hence  most  of  the  arithmetic  can  be 
followed  easily. 


6.  PROGRAMMING  THE  t2exp(-3t)  CASE 

A  simple  analysis  of  the  limiting  behavior  of  equation  (31)  shows 

that  as  r  — oo,  w  —  t2e"’  ,  ~  MAX — *  2’  and  WMAX  ‘  4^e  =0"54134 
When  5  -»  0,  w  —  0,  and  moreover,  from  (37),  when  5«  1, 


2*(sin?-  -  3?  cos  »*)+  ?a(b  +  4r  +  T2)e" 


(38) 


which  is  a  sinusoid  with  a  small  contribution  from  the  transient  term. 


It  is  clear  then  that  for  large  ?  we  may  be  certain  to  find  the 
maximax  upon  scanning  the  domain  of  §  between  1  and  3;  likewise,  for 
small  t,  we  may  expect  the  first  peak  to  be  the  maximax.  The 
situation  for  intermediate  values  of  §  is  not  clarified  by  a  brief 
additional  analysis:  e.g..  the  transient  term  has  no  relative  maximum 
at  the  level  of  approximation  given  in  equation  (38).  Practically,  this 
question  is  readily  solved  by  the  analog  computer.  The  results  of 
this  study  (details  are  furnished  in  section  6.  1)  show  that  for 
c  <  3,  the  first-peak  criterion  is  valid.  Also  the  criterion  for  f  >  3 
is  that  1  <  tmax  <  3. 

Based  on  these  two  criteria,  a  straightforward  digital  program 
was  written  in  FORTRAN  to  calculate  and  plot  the  shock  spectrum. 

A  short  explanation  of  the  program  is  given  in  section  6.  2. 


6.  1  Analog  Computer  Analysis 

An  elementary  circuit  is  sufficient  for  generating  w(t ).  All  the  inter 
mediate  functions  required  can  be  generated  by  linear  equipment 
(ref  9).  The  unsealed  equations  are: 
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Gj  =  e'T 

Gi'  =  -kGj 
G?  =  k(Gj  -  G?) 
Ga'  =  k(2Gp  -  G3) 
Hi  =  k§°(G3  -  w) 
w '  =  kHi 


k  =  dT / ds 


scaYed  cl«uiViriIvenainefimrer4SPfCl  1°  S'  the  math'"e  time.  The 
scaling  factors  that  depend  on  5  ’Ytabl  ^  ^  a"d  ^  are  unknown 
these  factors  is  given  in  appendix  C  Th'  g‘Vlng  8uitabll>  values  of 
a-  speed  of  the  recorde AVt^^  ™ the 


high  values  of  r  are  glveYiYfigYrY  5/^  i0W’  intermediate,  and 


these  response  curvf  AY^eYtYthe11  us^YY*  equalion  fading  to 
available  °n  the  EAl  231R  analo  r«-*petitive  operation 

used  here.  )  In  this  P-POP  technique^  (The  131R*as 

be  solved  in  fast  time,  and  the  parlnwirV Jlfferemial  equation  would 
dynamic  range  of  5  and  a  fortiori  *  L  ,  ‘  SWept  ln  slow  time.  The 
the  analog  computer;  scaling  problems  Unfor,unately  loo  great  for 
no  better  than  divide  the  t'-axfs  inlaT  ?”,e*  II  aPPears  one  can  do 
separately.  This  naturally  SCak‘  each  Action 

the  method,  which  was  therefore  not  pursued.  aUrac,ive"es*  °f 


— *-Z  DiSital  Computer  Program 


sented  in  figured.  bYaYYY*  th*  Sh°Ck  sPuctrum  is  pre¬ 

calculation  and  search  for  a  ms  °Ch  statement  after  DO  29  allows 
ing  the  maximaxes  inV&VT  «T<  ^  J  <  ’  <  ^  ?  >  3  «or- 
chosen,  which,  starting  from  -  -  0  ’  a  sumiglobal  algorithm  is 

m  fairiy  large  steps,  reverses  toosZ?**  f8t „the 
esweeps  and  repeats  until  the  interval  ’  CUtS  the  lnterval  by  five 

WMAV  thC  °ther  path  o{  the  Program  a  V*  r°Ughly  etiual  to  those’ 
WMAX.  The  SCAUT  subroutine  ^.’t.  Y  St0ring  ,he  results  in 

™surv  Pitting,  scales  TAU  and  WMAX  A’  Hau8ner’:’  for  simple 
tape  (in  low  density).  MAX,  and  writes  them  on 

,or  *u  ih*  <«». ...  . . . 


"  U°  "  ■CV"‘“  *Ub™"“~*  b»  HDL  , 


nternal  use. 
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7.  RESULTS 


The  chief  results  of  this  study  are  the  shock  spectra  shown  in  figures 
7  through  12.  As  can  be  seen,  figure  7  illustrates  a  family  of  forced- 
vibration  spectra  for  different  values  of  R  corresponding  to  the  forcing 
function  of  figure  la.  Figure  8  is  the  analogous  family  corresponding 
to  figure  lb.  Figures  9  and  10  are  the  free- vibration  spectra  for  the 
same  two  cases.  Figure  1  1  shows  a  lorced-  and  a  free- v, oration 
spectrum  on  tne  safne  graph  for  R  "  0.  01.  These  are  the  upper 
curves  in  figures  7  and  9.  Finally,  figure  12  is  the  spectrum  for  the 
forcing  function  of  figure  lc.  These  results  have  not  been  observed 
in  the  literature. 

The  important  features  of  figures  7  through  11  are  the  resonance 
peaks  near  y  =  1.5,  10,  20,  and  30.  If  we  were  to  "continue"  the 
half-sine  pulse  to  a  continuous  monotonic  frequency,  then  y  would 
measure  the  natural  frequency  of  the  system  with  respect  to  the 
former  frequency.  When  there  are  five  positive  impulses  per  half¬ 
sine  pulse,  the  system  responds  to  this  repetition  rate  of  ten  per 
"cycle"  in  giving  the  resonances  at  10,  20,  and  30  with  negligible 
augmentation  of  the  major  peak  of  the  pure  half-sine  pulse  spectrum 
(R  =  0).  Additional  energy  is  fed  into  the  system  when  interspersed, 
negative  impulses  are  added;  the  amplitude  of  the  spectrum  is  thus 
larger  except  for  suppression  at  values  of  y  at  even  multiples  of  10. 
The  suppression  is  due  to  the  cancellation  of  the  kinetic  energy  of 
the  system  oscillations  by  the  forcing  function.  This  happens  when 
the  system  has  just  completed  an  oscillation  cycle  and  is  then 
"s.ruck"  by  a  negative  impulse.  The  velocity  at  this  moment  is 
maximum  in  the  positive  direction,  but  the  force  is  in  the  negative 
direction. 

In  figure  11,  we  see  how  the  residual  spectrum  may  predominate 
for  certain  values  of  y,  in  this  case  for  y  <  1;  one  must  in  general 
then  consider  both  the  forced  and  the  free  era. 

The  tabular  data,  not  given  here,  show  that  in  no  case  is  the 
negative  maximax  larger  than  the  positive  for  the  forced  vibrations. 
(They  are  equal  in  the  residual  era.  )  The  trend  of  the  data,  however, 
strongly  suggests  that  for  resonances  at  larger  y,  the  negative 
maximaxes  may  become  larger  than  the  positive.  Moreover,  some 
results  generated  for  R  =  0.  2  with  the  earlier  program  showed  that 
for  the  alternating  impulse  case,  negative  maximaxes  were  larger 
near  and  at  the  resonances  for  y  =  10  and  30.  Thus,  although  not 
investigated  here,  negative  maximaxes  may  be  very  significant 
for  large  y  and  large  R. 
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DIMENSION  TAUIIOOOI.TAUSI 1000 1 , SXNTAUI lOOOf , *1 1 JSO) , WMAX 1 350 ) , 

1  AI3  JOI.AU350I.A2I  3501,811  3301(821130(830501 

1  FORMAT 1 141 5) 

11  FCRMATIBFIO.SI 

2  FORMAT ( IH1 1 

12  FCAMAT I  OCX?  I  TX2HX 1 13X4FWMAX4X 1 1 // I 
22  FORMAT!  I  2CX2IOPF  10.2,  1P220.7))) 

32  FCftMATimSftll  •  I4,?X5FJJ  ■  131 

42  FORMAT 13TX6HCX l  •  F5.2,SXTM0TAU  •  FT.4/1TX4NDELTAM  •  F5.2.5X4H2  > 
l  FT. 4///) 

READ  1 5, 1  I  It  tJJ. I  If I?iNXfNYt I DENT, ICOCE 
ME AO  15 ,1 1 1  0Xl,DTAU,0ELTAM,2 
C 

C  II  •  NUNOEA  CF  VALUES  CONFUTEO  )F  t ACM  RESPONSE  FUNCTION 

C  JJ  •  NUMBER  CF  POINTS  IN  THE  SPECTALM 

C  11  •  PLOTTING  OPTICN  (PLOT  WHEN  1 1 - I ) 

C  12  •  WAITING  CPTION  IWR1TE  TABLE  OF  RESULTS  WHEN  12*11 

C  0X1  <  XI  INCREMENT 

C  DTAU  •  TAL  INCREMENT 

C  DELTAM  •  INITIAL  INCREMENT  IN  TAM  REPRESENTING  TAU 

C  DELTAN  •  LPOATEO  INCREMENT 

C  l  *  eRROA  CRITERION  IN  TAJ 

C 

C  OC  L UOP  9  GENERATES  ARRAVS  BASSO  ON  I .LT.TAU.LT.? 

C 

DC  9  1*1.11 

FI  •  I 

TAUI I)  •  F|*CTAU  4  1. 

TAUSIII  •  TAL(I)fTAUII) 

EXMTAUII)  •  EXPI-TAUI 1 1  I 
4  CCNTINUC 

OC  19  J*1,JJ 
WFAXIJ)  •  0..0 
FJ  .  J 

XI(J)  *  FJ.OXI 
X1S  •  XIIJMXIU) 

ONLPKS  •  1.  ♦  X*S 
ONPXSS  •  CNEPXS*LNEPXS 
ONPXSC  *  CNPXSS*CNEPXS 
A1JI  •  2.«XI(  JI/ONPXSC 
AlIJ)  •  1.  -  5.4XIS 
A2IJ)  •  -«l(  JKO.-X1SI 
B1IJI  •  •*U)*A?(J) 

82IJI  •  4,*XIS/0NPXSS 
B3IJI  •  XIS/CNEPXS 
19  CONTINUE 

OC  29  J.I.JJ 

IFIXIUl.LT. 3.01  GO  TO  100 
DC  J'»  1*|.II 

SIMIT  ■  S  INI  X  I  U I  *TAOI  1)1 
CCSXIT  •  CCS  I  X 1 1 J )*TAUI 1 1  ) 

W  •  AU)*(AI(  JKSINXIT  «  A2U)*C0SXiri  4  ID1IJ)  4  R2UI*TAUIII  4 
l  F3U)*TAUSIIII*CXMTAUI  I) 

IF!  W.GE.kPAXI J) )  WKAYIJI** 

39  CONTINUE 
GO  T  Ci  2< 

100  T**>  •  O.u 

DELTAS  *  CELTAM 
WFULOl  *  C.O 
80  VMULI'2  «  C.O 

Figure  6.  Program  TSEXP. 
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60  TAN  •  TAM  ♦  CELT AN 
XITAM  .  XIIJl.TAh 

**«■•»•“»«*"«»»♦  l.u, 

IFIW.LT.MFCLCI)  CO  TO  50 
NM0L02  •  kHOLCl 
MH0L01  •  k 
GC  TO  60 

so  ifidlltan.lt. ;>  r.o  to  to 

TAM  •  TAM  -  2.«0ELTAN 
MMOL  01  •  kM0LC2 
CELT  AN  «  CSLTAN/5. 

GC  TO  60 

70  MMAX(J)  •  MHCL01 
24  CONTINUE 

lFIIl.NE.il  GC  TC  10 

*CU,T,XI  «MMAX,JJ,NX.NY,  I02NT,  ICCOt  I 
“EWI  >10  9 

10  IFII2.NE.il  STOP 
JCOUNT  •  C 
MR  IT  1 1 6. 2 1 
MR !T£ I  to. 32  I  II, jj 
WRITE  I  to, 62  I  CXI.CTAU.OELTAM.I 
JO  J1  ■  1  ♦  I00» JCTUIT 
J2  ■  J1  ♦  A9 
WRITEI6, 12  ) 

JCWNT4:2jcnuNTli,i""**<J,,*MJ*!0,,N''AX,J*5°,,J*J,,-m 

I FIJJ.L 3. 100* JCOUNT I  GO  TO  20 

WRITEI6.2) 

GC  Til  30 
20  STOP 
EkC 


B2(  JUTAM 


Figure  6.  Program  TSEXP. 
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Figure  7.  Forced  vibration  spectra  of  five  positive  impulses  on  a  half¬ 
sine  pulse  driving  an  ODOF  system. 
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,  Free  vibration  spectra  of  five  positive  impulses  on  a  half¬ 
sine  pulse  driving  an  ODOF  system. 


Figure  9 
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- 

Figure  11 


>  *  I  j  *  ID  IB  14 

.  Comparison  of  forced  and  free  vibration  spectra  for  ten 
alternating  impulses;  R  =  0.01. 
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APPENDIX  A.  —  EVALUATION  OF  COEFFICIENTS  OF  PARTIAL 

FRACTION  fcXf>ASi316Fi — — — - 


If  in  (33)  we  multiply  by  s  and  let  s  — oo, 
A  +  C  =  0 
If  we  let  s  —  0, 

B/f  +  C  +  D  +  E  =  \/f 
If  we  multiply  (33)  by  (s  4  1  )3  and  let  s  —  -1, 

E  =  1/(1  +  f) 

If  we  multiply  by  s'  +  and  let  s  —  ic , 

Ai  r  +  B  =  1/(1  +  iS)3 
Equation  (A-4)  yields 

(Aif  +  B)(l  +  3i*  -  if  -  e?a)  *  1 


(A-  1 ) 
(A-2) 
(A-3) 
(A-4) 


(A-5) 


which  upon  expansion  and  equating  real  and  imaginary  parts  yields  the 
pair  of  equations 


and 


( 1  -  3^'  )B  -  f(3  -  f)A  =  1 
(3  -  f  )B  +  (1  -  3f)A  =  0 


(A-6) 

(A-7) 


The  determinant  of  the  coefficient  matrix  reduces  to  ( 1  +  *L)3,  The 
solutions  are  readily  found  to  be  those  given  in  (34a)  and(34b).  Finally 
D  is  obtained  from  (A-2) 
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One  defines  T 

I  =J  v2e'vsin5(T-v)dv 

o 


J  =  v?e"v  cos  f  (t -v)  dv 

O 

T 

K=  j  ve“vsin  ?(t-v)  dv 


L,=J”  ve‘vcos'(*-v)dv 


M=  j  e"vsin  £(t-v)  dv 


(B-l) 
(B-2) 
(B-3) 
( B  —  4 ) 
( B-  5) 


N  e‘vcos  f('-v)dv  (B-6) 

where  w  =  '  I,  Integration  by  parts,  retaining  e*v  under  the  differential 
each  time,  yields  the  relations 

I  =  2K  -  *J  (B-7) 

J  =  -  ' e  +  2L  +  ll  (B-8) 

K  =  M  -  FL  (B-9) 

L-  =  -*e“~  +  N  +  ?K  (B-10) 

M  =  sin  -  'N  (B-  1 1) 

N  =  -e"'  +  cos'-  +  (B-12) 

whose  solution  for  I  may  be  given  as 
1(1  + O'"  =  :(l  +  ?:T':e-'  +  4?(l  +  f>e-~  +  2'(3-'  )e‘*  + 

2(l-3Fr)sin;  -  2  '(  3  -  )cos  '  •  (B-l  3 ) 
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APPENDIX  C. 


SCALE  FACTORS  OF  ANALOG  CIRCUIT 


Table  Cl  summarizes  the  essential,  approximate  scaling  data 
obtained  by  trial  using  the  circuit  in  figure  4.  The  potentiometer  num¬ 
bers  in  the  five  columns  at  the  right  are,  of  course,  arbitrary.  The 
circled  numbers  denote  gains  on  the  EAI  131-R  analog  computer.  The 
values  and  gains  are  easily  translated  into  corresponding  values  for 
computers  having  different  gains  available. 


Table  Cl.  Scale  Factors  and  Potentiometer  Settings  for  Analog  Circuit 


? 

k' 

■  CW 

Si 

42 

44 

ft - 

48 

52 

0.2 

2 

500 

2500 

.5437  j5 

.  1259  0 

.  5413  © 

.4 

.  4 

0.5 

1 

125 

250 

.2718  l7j 

.6796  2' 

.  3383 

.  5 

.  5 

1.0 

1 

100 

100 

.2718  fb 

.6796  (2 

.5413 

T 

1 

n/3 

1 

100 

100 

.2718  (TO 

. b796  © 

.8120 :2 

(T) 

.6© 

2.0 

1 

100 

100 

.2718  ft 

.6796  (£' 

.4331  © 

ft 

.8  'S 

3.0 

1 

100 

100 

.2718  Q 

.6796  -J 

.9744© 

X 

.9  fb 

4.0 

1 

100 

100 

.2718  0 

.6796  T 

.8661  fd 

X 

5.0 

0.2 

125 

100 

.5437 

.2718 

.5413  © 

.25 

,8  © 

10.0 

0.  1 

125 

100 

.2718 

.1359 

.5413  0 

.  125 

.8  Jj) 
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APPENDIX  D. —  TABLE  OF  SYMBOLS 


m 

k 

x 

t 

F(t) 

Fo 

e 


mass 

spring  constant  (on  analog  diagram,  used  as  time  scale 
factor) 

displacement 

time 

forcing  function 

maximum  value  of  half-sine  pulse  forcing  function 
(excluding  impulses) 

given  decay  constant  of  one  forcing  function 
natural  frequency  of  system 
static  displacement 


U  xs 
lc 

t 


lc 


'  F 


TM 

UM 

F; 


4 

ill  -  — —  u 

e' 


e  =  Vrr 


response  factor 

duration  of  half-sine  pulse 
normalized  time 

angular  frequency  characteristic  of  forcing  function 

normalized  natural  frequency 

time  at  which  u  reaches  a  global  maximum 
value  of  u  at  its  global  maximum 
strength  of  Delta  impulse  function 

ratio  of  impulse  strength  to  half-sine  pulse  peak 

scaled  response  factor 

scaled  normalized  natural  frequency 


36 


unclassified 


DOCUMINT  CONTROL  DATA  .RAD 

■  •  OW'IIW. Ttw  activity  (CupttmH  1  1  ;aS .. aSSaelili 

*«*0*T  SCCURITV  Cl  AMI  I 

Harry  Diamond  Laboratories  Unclassified 

Washington,  D.  C.,  20438  '*  •*0*" 


•  ■  NINORT  TITL« 


HOCK  SPECTRA 


4.  oi(cni4Tivt  NOTMfiy^afn^ana, ShUMnem, 

»U  r»o«m  ,flm  am,,  alii.  UliUl.  It;  SSS! - 

Paul  G.  Hershall 


nr  Oats 

January  1968 


a  nnojicr  no. 


..  DA-1N22Z901A301 
HDL  Proj:  39800 


Distribution  of  this  document  is  unlimited. 


UNCLASSIFIED 


UNCLASSIFIED 


Dirac  impulses 
Half-sine  pulse 

One-degree-of-freedom  system 
Idealized  forcing  functions 
Digital  programs 
Analog  computer  as  auxiliary 


